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The interaction of spatially small dipoles with nucleons, nuclei is calculated in the DGLAP approximation at 
the top of HERA energies and found to be close to the S'-channel unitarity limit in the case of the color octet 
dipoles. The DGLAP analyses of the current diffractive data appear to support this conclusion as they indicate 
a ~ 30 — 40% probability of the gluon induced diffraction for Q 2 ~ 4GeV 2 . The need for the high- precision 
measurements of the i-dependence of inclusive and exclusive diffraction for pinpointing higher twist effects in the 
gluon sector is emphasized. The eA collisions at HERA would provide a strong amplification of the gluon densities 
allowing to reach deep into the regime of nonlinear QCD evolution. Connection between the leading twist nuclear 
shadowing and leading twist diffraction in ep scattering is explained. The presented model independent results for 
the nuclear shadowing for light nuclei indicate much larger shadowing for the gluon sector than for the sea quark 
sector. It is argued that HERA in eA mode would be able to discover a number of new phenomena including 
large gluon shadowing, large nonlinearities in parton evolution, small x color transparency in the vector meson 
production followed by color opacity at x < .005, large probability of inclusive diffraction. Implications for the 
nucleus-nucleus collisions at LHC are discussed as well. 



1. INTRODUCTION 

Fascination with small x - physics can be traced 
back to the old questions - why is the cross section 
of nucleon-nucleon interactions at high energies 
large and slowly growing with si Can the interac- 
tion of a heavy onium state with hadrons remain 
weak even at very high energies? Is it possible to 
exploit the increase of the parton distributions at 
small x to produce superdense hadron matter in 
the laboratory? 

A unique advantage of DIS for studying dy- 
namics of strong interactions is the possibility to 
tune the resolution scale. Nuclei add another di- 
mension - possibility to vary the thickness of the 
target. Hence it is possible to investigate the 
transition from soft nonperturbative dynamics to 
perturbative dynamics: up to what virtualities 
(or down to which distances) is nonperturbative 
QCD important? In addition, how does the tran- 
sition depend on the incident energy? 
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The fundamental questions one can address in 
DIS at small x are: How large is the cross section 
of interaction of small "hadrons" (oniums) with 
nucleons /nuclei at high energies? Is a kind of the 
Froissart bound for DIS close to the HERA kine- 
matics for the gluon induced processes?, Would 
nuclei be transparent to a high-energy J/ip,Y? 
How far down in x it is necessary to go for the 
observation of a break down of the linear DGLAP 
evolution equations? What are implications of 
nonlinear QCD evolution for the physics of AA 
collisions? 

In this talk we first focus on ep scattering. We 
discuss the interaction of spatially small color sin- 
glet clusters with hadrons at high energies. The 
S'-channel unitarity constraint on this interaction 
implies that, in the case of color octet dipoles 
interaction may already be close to the unitar- 
ity limit already at higher HERA energies, im- 
plying a large probability of diffraction when the 
hard process is induced by the interaction with a 
gluon. The use of the QCD factorization theorem 
[ |l| for inclusive diffraction allows the extraction 
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of this probability from the current analyses of 
the HERA data. The probability is indeed large 
~ 30 — 40% in a wide range of Q 2 which is much 
larger than the measured probability of the quark 
induced diffraction. We found that breakdown of 
the DGLAP approximation for the gluon chan- 
nel (higher twist effects becoming comparable to 
the leading twist effects) is related to large parton 
densities at low Q scale rather than to the BFKL 
type dynamics [ ^) . It should be accompanied by 
large fluctuations of gluon densities. Further di- 
rections of studying the gluon dynamics at small 
x at HERA are outlined. 

In the second part of the talk, we explain how 
the study of electron-nucleus collisions at HERA 
can add to the understanding small x dynamics. 
We emphasize that a number of expected phe- 
nomena follow from the basic features of the small 
x dynamics. In particular, the application of the 
S-channel unitarity shows that nonlinear effects 
in the gluon channel already are large. The Gri- 
bov black body limit for the cross section of DIS 
[ || is found to be a good guess for gluon inter- 
actions in the HERA kinematics. Combining the 
Gribov theory of shadowing and the QCD factor- 
ization theorem for diffraction allows the deriva- 
tion of the model independent expressions for the 
leading twist quark and gluon shadowing for the 
scattering off light nuclei. An enhancement of 
gluon shadowing, compared to quark shadowing, 
is predicted. It is concluded that nuclear matter 
at HERA energies should be opaque to the small 
size dipoles, leading to break down of color trans- 
parency regime at small x (x < 0.005 for J /if) 
photoproduction). Large gluon shadowing would 
lead to significant changes of the first stage of the 
heavy ion collisions at LHC. Really account of nu- 
clear shadowing within the conventional convolu- 
tion formulae leads to a strong (a factor ~ 10) 
suppression of the yields of minijets. 

2. ELECTRON-PROTON PHYSICS 

2.1. The interaction of small color singlets 
with hadrons 

At small x the average longitudinal distances 
(Ioffc distances) in the correlator of the e.m. 
currents which determine F%p{x,Q ) are l co h ~ 



— - — . At HERA l rn h can reach 10 3 fm. Hence in 
the target rest frame the small x processes occur 
in three steps. Initially the 7* transforms into 
pair. Next, the 33 color dipole of transverse 
size b interacts with the target. This cross section 
rapidly grows with increase of energy [ ||, ^] : 

2 \ 
°%%{E mc ) = ^-b 2 a s (Q 2 )xG N (x,Q 2 ~ -), (1) 

where \{x « lO" 3 ) « 9,x = 2m %_ ■ 

For color octet dipole the cross section is en- 
hanced by the ratio of Casimir operators of color 
group SU(3) C - C F (8)/C F (3) = 9/4 [ |,|: 

"TfAEinc) = ^b 2 a s (Q 2 )xG N (x,Q 2 ~ A), ( 2 ) 

These equations are another form of the QCD 
evolution equations. For the fit to the energy 
dependence ^| a(s,Q 2 ) cx s( a (Q one finds 

that due to increase of the steepness of the gluon 
density with virtuality, a{Q 2 ) gradually increases 
with Q 2 : 

a{AGeV 2 ) w 0.2, a(A0GeV 2 ) w 0.4. (3) 

2.2. Pomeron: one, two, ... too many? 

Satan: "You were selling spoiled sturgeon fish. " 
Merchant: "No - sturgeon was of regular quality" 
Satan: "Sturgeon could only be the prime quality 
there is no other one. " 

"Master and Margaret" M.Bulgakov 
Eqs.[l],|],|] clearly demonstrate a qualitative dif- 
ference of the strong interaction dynamics in 
DIS at small x from the soft physics of the 
vacuum exchange. This object was derived by 
V. Gribov within the high-energy Regge theory as 
the rightmost singularity in the complex angu- 
lar momentum plane. It was named Pomeron 
by M.Gell-Mann after I.Pomeranchuk who was 
the first to predict the total cross sections of 
particle (antiparticle)- hadron scattering to be 
equal at asymptotic energies. The key feature 
of the Regge trajectories is their universality - 
independence of the trajectories on the process. 

2 Such a fit is useful in the practical applications although 
QCD predicts a behavior <x a + felni + cln 2 A for the 
HERA energy range since radiation of oc 1-2 hard gluons 
is possible in the multiRegge kinematics at HERA. 
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Clearly the behavior of the cross sections, in the 
kinematics where DGLAP describes the Q 2 evo- 
lution of the parton densities, violates such uni- 
versality in a gross way. 

It is sometimes suggested that one can save the 
idea of the Pomeron in DIS processes by introduc- 
ing two Pomeron trajectories with different inter- 
cepts. Clearly a smooth change of the effective 
power with Q 2 indicates that this hypothesis is 
not equivalent to eqjl]. 

Another manifestation of the breakdown of the 
universality is the change of the pattern of the 
i-dependence of the quasi two-body processes ex- 
pected in pQCD [ [jj. Let us consider for example 
the process 7* + p — > V + p where V is a vector 
meson. 

In the soft regime the t dependence of the dif- 
ferential cross section changes with s as 

/(s,0)cK*)exp(2a'*m(s/so). (4) 

In the impact parameter space this amounts to 
an increase of the radius of the projectile hadron 
with increase of the incident energy. This can be 
understood as due to multiperipheral structure 
of the soft process - fast partons of the projectile 
decay into slower partons with a random change 
of the position in the impact parameter space of 
b ~ -jr- where kt is the soft scale (transverse mo- 
mentum of the emitted parton). The resulting 
random walk in the impact parameter space leads 
to a growth of (6 2 ) proportional to the number 
of steps, that is to the length of the rapidity in- 
terval. Hence the size grows oc Ay, leading to the 
shrinkage of the cone given by eq.|| [ || . 

In the hard regime the parton (gluon) emissions 
are governed by the hard scale (Fig.l). Hence the 
random walk produces much smaller increase of 
(b 2 ): a ha rd < a so ft- 

Therefore with an increase of Q 2 one expects 
three phenomena to happen [ |t| : (i) a decrease of 
the contribution of the upper end of the ladder 
to the slope, B (f(t) = exp(Bt)) which should 
tend to the limiting value determined by the two- 
gluon form factor of the nucleon - experimen- 
tally £?2 9 ~ 4 GeV~ 2 . (ii) Gradual disappear- 




Figure 1. Electroproduction of J ftp. The rapidity 
interval occupied by hard exchanges 




In s 



Figure 2. Sketch of the Q 2 and energy depen- 
dence of the slope of the vector meson production 



ance of the energy dependence of the slope, (hi) 
A fast increase of the cross section with energy 
~ |a;G(a;, Q 2 ) | oc ^4. However this pattern can 
be sustained only until the transverse momenta 
in the perturbative ladder degrade to the level 
corresponding to the soft physics. At this point 
the further steps in the ladder will be soft and 
will lead to increase of the slope (see the sketch 
Fig. 2) and lead to the values of a' which are likely 
to approach, at very high energies, values similar 
to those for the soft physics. 

Obviously the discussed pattern is qualitatively 
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different from the soft Pomeron type physics. 
Hence, one should be extremely careful with the 
use of the word "Pomeron" in the DIS regime. 

Current HERA data confirm the prediction of 
the convergence of the slope of p-meson produc- 
tion at large Q 2 to that for J/ip (section 2.5). 
However the HERA data are not sufficiently ac- 
curate so far to determine a' for DIS processes 
in the HERA energy range and the change of the 
slope between the HERA range and the fixed tar- 
get experiment range is still a matter of debates 
[ |9|, [IFJ. Nevertheless one can obtain clear in- 
dications of nonuniversality of a' by combining 
HERA and SLAC linac data. Indeed, assuming 
for a moment that a' is universal we would find 
that B«4 GeV" 2 for W = 200 GeV as measured 
at HERA would lead to B « 1 GeV" 2 at W ~ 6 
GeV, while the SLAC linac data at that energy 
correspond to B ss 3 GeV~ 2 . 

Hence we conclude that an accurate mea- 
surement of the energy dependence of the t- 
dependence would provide a unique way to de- 
termine the transition from soft to hard physics 
and provide a direct evidence for non-Regge type 
behavior of the interaction in the vacuum channel 
in DIS. 

2.3. Can cross sections grow rapidly for- 
ever? 

A fast increase of the inelastic cross section of 
the small dipole -nucleon interaction cannot con- 
tinue indefinitely. The black body limit: (Jdiff < 
<j tot 1 1 indicates that this is impossible. 

Nucleon case 

The analysis is much simpler than in the case 
of soft interactions since the radius of the system 
practically does not grow with energy (see discus- 
sion in the previous section and Fig. 6). 

It follows from the optic theorem that [ |ll[ : 
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(5) 



where r\ is the ratio of real and imaginary parts of 
the amplitude. Using the Gribov-Migdal relation 
for rj: 



_ 7T d In a tot 
V ~ 2 dins 



and T} > 0.4 for Q 2 ~ 10 GeV 2 from eqsjl], | we 
find a kind of the Froissart bound: 



o-tot < — — - ~ 35m6, 



1 



(7) 



for Q 2 = 10 GeV 2 and therefore <x, ne z < 18mb. 
The limit becomes tighter with increase of Q 2 due 
to increase of rj with Q 2 . 

To illustrate numerically to what extent this 
limit is relevant for dynamics at HERA and be- 
yond let us consider Q 2 ~ 10 GeV 2 and x ~ 10~ 4 
and take xGm ~ 20 and as = -25. Clearly such 
an estimate carries substantial uncertainties, in 
particular because eqs.(|lj|^) were derived in the 
leading a s log Q 2 approximation. The cross sec- 
tion of interaction of the color triplet dipoles with 
nucleon is still much smaller than the discussed 
limit: Umei ~ 7 mb. Using cq.(^) we also get 
a " 1 w 0.1. At the same time for the color 

find 

octet dipole eq.(gj) gives <Jinei(x — 10 _4 ,Q 2 = 
10GeV 2 ) w 14 mb. This is pretty close to the 
limit and corresponds to: a e i/a to t ~ 1/3. This 
indicates that decomposition of cross section over 
twists in this situation becomes unreliable, also 
one should expect a large probability of diffrac- 
tion in the gluon induced processes. Applying the 
AGK rules [ [T^] for the diagrams corresponding 
to the higher twist terms (like attachment of four 
gluons to the dipole) we find that they lead to 
final states with double, triple, ... multiplicity. 
This provides an evidence for the large fluctua- 
tions of the parton densities in the regime close 
to the unitarity limit f] 

It is worth emphasizing that a large value of 
xGn(x, Q 2 ), which results in a breakdown of the 
linear evolution equations, was generated within 
DGLAP predominantly due to logQ 2 evolution 
rather than solely due to In 1 jx effects which 
would be the BFKL approximation. Qualita- 
tively the large gluon density is due to the pres- 
ence of large number of gluons already on a low 
resolution scale (gluon carry about the same mo- 



(G) 



3 It is worth noting that there is no direct relation between 
the higher twist effects in diffraction and in the total cross 
of DIS. In the diffraction case we could use the optic the- 
orem to single out the leading term in the scattering of 
small clusters. In the case of the total cross section sev- 
eral competing terms enter with opposite signs, cf. [ Q . 
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mentum as quarks at this scale) and they gener- 
ate a lot of gluons in the evolution. (At HERA 
kinematics xG(x,Q 2 ) « 20 but only 1-3 gluons 
are due to multiRegge kinematics). When ener- 
gies approach the unitarity limit the increase of 
<t(&, x) should slow down to a rate which may be 
rather similar to the case of soft physics. One 
may expect that in this limit a soft regime will 
be reached with a genuine universal Pomeron. 

Nucleus case. In the case of a sufficiently heavy 
nucleus the inelastic cross section cannot exceed 
the geometric size of the nucleus: 

<-a(M) = ^b 2 a s xG A {x,Q 2 = A/6 2 ) < nR 2 A (8) 

Hence xG A (x, Q 2 = 10)/ A < 120A- 1 / 3 . For cen- 
tral impact parameters the inequality is a factor 
of 1.5 times stronger: 

xG cent.imp.par.^ Q 2 = 1Q y A < 80A -l/3_ (Q) 

For the interactions induced by the color octet 
dipole the constraint is a factor of 9/4 stronger: 

xG c A entr imp par (x,Q 2 = 10) /A < 35A- 1 / 3 , (10) 

which for A = 200 is substantially smaller than 
xGat(10- 4 ,10) - 20. 

Thus we conclude that at the HERA kinemat- 
ics quark induced interactions do not reach the 
unitarity limit but perhaps the gluon induced in- 
teractions are close to it (estimates at lower Q 
are more uncertain but the trend is that the ap- 
proach to the unitarity limit should occur earlier 
for smaller Q). If the gluon shadowing were small, 
the quark induced interaction at small impact pa- 
rameters for heavy nuclei would be close to the 
unitarity limit for Q 2 ~ 10 GeV 2 leading to large 
nonlinear effects in the QCD evolution for quark 
sector. 

In the case of heavy nuclei the unitarity con- 
straint definitely points to a modification of the 
dynamics for the gluon induced interaction even 
if the gluon shadowing is large. 

Extension of the x-range by two orders of mag- 
nitude at TESLA-HERA collider would corre- 
spond to an increase of the gluon densities by 
a factor of 3 for Q 2 = 10 GeV 2 . It will defi- 
nitely bring quark interactions at this scale into 



the region where DGLAP will break down. For 
the gluon-induced interactions it would allow the 
exploration of a non-DGLAP hard dynamics over 
two orders of magnitude in x in the kinematics 
where a s is small while the fluctuations of parton 
densities are large. 

2.4. Inclusive diffraction 

We discussed above that a large probability of 
diffraction in DIS would be a clear signal for the 
onset of a new nonlinear dynamics. Since this 
characteristic of DIS can be directly measured, 
it provides a more direct probe than analysis of 
inclusive structure functions which depend on the 
nonperturbative boundary conditions which are 
not constrained by the theory strongly enough. 

The current situation with diffraction at HERA 
is already very intriguing. One expects, based 
on the QCD factorization theorem for diffractive 
processes [ |) , that the diffractive structure func- 
tions fj-'(P,Q 2 ,xp,t) satisfy DGLAP evolution 
equations. Here [3 — x/xp. The data seem to 
be consistent with an early scaling for Q 2 > 4 
GeV 2 for most of j3 range and a large probability 
of diffraction - P q ~ 10% for j*p scattering poses 
rather strong constrains on the dynamics. 

This experimental observation is highly non 
trivial. Indeed, if the DGLAP evolution for the 
fragmentation would be valid for Q > Qa = 
0.7 GeV where quark and gluon parton densi- 
ties tend to zero at x — > (the GRV scenario) 
the rapidity gap probability in DIS would be 
« i 5 , k 0.1 observed at HERA. (Naively, if one 
assumes factorization for the inclusive cross sec- 
tions, the assumption of an early factorization for 
the leading baryon production should be reason- 
able as well. Indeed, the overlapping integral for 
the leading partons (present at the scale Qq ) with 
the emitted partons to form a low p t leading nu- 
cleon should be very small since the partons which 
are generated in the evolution have p t > Qo and 
rather small light-cone fractions ). 

Another interesting feature of the current data 
(shared by the diffractive data from the proton 
colliders) is a very important role of the gluons in 
the diffractive events. In the parlor of diffractive 
community - the "Pomeron" is predominantly 
build of gluons. To quantify this statement it is 
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Figure 3. Probabilities of rapidity gap events for 
quark and gluon induced diffraction calculated 
using the best fit of [ [f4| . 



convenient to define the probability of diffraction 
for the action of the hard probe which couples to 
a parton j: 



Pj(x,Q 2 ) 



I f^(^,Q 2 ,x ]P ,t)dtdxj P 
fj(x,Q 2 ) 
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If the interaction in the gluon sector at small 
x reaches strengths close to the unitarity limit 
we should expect that P„ is rather close to 



expect unc^ ± g 

1/2 and much larger than 

for rr inel 

color octet dipole—N 

/O 2 — in n^XT? 



Using eq.[ 



we calculate P g (x 



Q 2 = 10 GeV 2 ) « 1/3. Indeed, using the 
global fit of Alvero, Collins and Whitmore [ [l4| 
we find (Fig.3): P g > P q , and P g {x < 10~ 3 ) « 
0.4(0.3) for Q 2 = 4(10) GeV 2 ! 

ft is also interesting that diffraction when pa- 
rameterized in terms of the effective Pomeron tra- 
jectory corresponds to ap(t — 0) = 1.15 — 1.2 
which is somewhat larger than what is usually 
observed in the soft domain. 

A word of caution: the direct measurements of 
fg° are available for a limited range of xp and 
for rather large Q 2 . "Backward" evolution ne- 
glects higher twist effects which may slow down 



this evolution and reduce P g at Q ~ Qo- How- 
ever backward evolution should be a small effect 



for Q 2 



15 GeV 2 where P„ ~ 0.3. So we con- 



clude that the data confirm the expectation that 
interactions in the gluon sector are strong. They 
are not corresponding to the black limit since P g 
is substantially smaller than 0.5. Also, the mea- 
sured slope of diffraction suggests that configura- 
tions involved in the total cross section of diffrac- 
tion have a rather large size while the slope of the 
cross section of diffraction in the gluon induced 
processes was not measured so far. 

Due to the large distances involved in the 
diffraction we can think of diffraction as "elastic" 
scattering of qq or gg configurations of the tar- 
get (if we define the process at a given Q scale, 
changing the scale leads to a mixing of these con- 
figuration, and the addition of new intermediate 
states like qqg). Hence it is instructive to treat 
diffraction in a complementary S'-channel picture 
-the eigen states of the scattering matrix [ [lj|. 
Use of the optic theorem leads to 



a eff (x,Q ) 



.,. _ I6nda dlf M| t=0 



(12) 



and allows to extract the average cross section for 
the configurations which contribute to quark and 
gluon induced diffraction 



gluons ( 

eff 



(x= 10- 4 ,Q = 4) 



value of a 



gluons 
eff 



(x = 10- 



aff}{x,Q). We find 
30m6 and even larger 
2 = 2)- bOmb which 



indicates presence of superstrong interactions in 
the gluon sector (Fig.^). These strong interac- 
tions should be due to some hard dynamics of 
strong gluon fields since the 7+p — * J/ip+p data 
suggest that pQCD for these Q occupies most of 
the rapidity interval. 

This argument for superstrong gluon interac- 
tions at small x is complementary to the unitarity 
argument given above. 

Several avenues of studies are possible in ep 
scattering to investigate relative role of soft and 
hard physics in the gluon induced diffraction: de- 
tailed investigation of charm diffraction, diffrac- 
tive dijet production, energy dependence of the 
t-slope. The study of scattering off nuclei - nu- 
clear shadowing, diffraction off nuclei would be 
able to provide direct answers to these questions. 




(suppressing the z-dependences) 



Figure 4. The ir-dependence of <J e ff for quark 
and gluon induced diffraction. 



To summarize, studies of diffraction at HERA 
appear to confirm the conclusion of section 2.3 
that large nonlinear effects could be taking place 
in ep scattering for Q 2 < 4 — 10 GcV 2 in the gluon 
sector - the main problem is the lack of effective 
"gluonometers" (especially in diffraction) at this 
resolution scale. 

2.5. Vector meson production - a quest for 
a low Q gluonometer 

Over the last few years a lot of experimental 
and theoretical studies of the exclusive vector me- 
son production in D1S were performed. A proof of 
the QCD factorization theorem for exclusive me- 
son production by longitudinally polarized pho- 
tons [ |l6| has demonstrated that for large Q these 
processes may be calculated as rigorously as the 
leading twist processes. In the small x limit the 
amplitude of the process - A^l + p — ► V + p) 
can be written as a convolution of the light-cone 



wave function of the photon 



♦!<?<?) 



the scat- 



tering amplitude of the hadron state, a(qq~T), 
and the wave function of the vector meson, ipy. 
tr(&, s) is expressed through the skewed gluon den- 
sity of the target In impact parameter space 



4 At small x and large Q 2 the skewed parton distributions 
are calculable in QCD as a solution of the QCD evolution 
equation with the skewed kernel and the conventional par- 
ton distribution at low scale as the nonperturbative input. 



A 



d 2 bip~/*,L(b)cr(b, s)ipy(b). 



(13) 



The parameter free leading twist answer for the 
process is 



r1n L 



l(Q)vv 



dt 

2 \f! 



t=0 = 127T Ty^ e + e - My 



'2 din x 



xG skewed, 



a E AiQ 6 N 2 



(14) 



Here, T v ^ e + 

f dz d 

1 J ~ 



is the decay width of V — > e + e ; 

J>v(z,k t ) , „ 2 

„ T — — ; > 3 lor Q z — > 00. 

2 J dzcPkt *v(»,*«) 

In the leading twist b=0 in ipy(b). Finite b ef- 
fects in the meson wave function in the overlap in- 



tegral fll3[ ) appear to be one of the major sources 
of higher twist effects at small x. They strongly 
reduce the cross section at intermediate Q 2 . How- 
ever the transverse sizes essential in the ampli- 
tude may remain small, down to rather small Q 2 
(FigJ|). Hence it is likely that the study of p, (p- 
meson production at intermediate Q 2 ~ 5 — 10 
GeV 2 will allow the study of the gluon densities 
at Q 2 down to 2-3 GeV 2 . The best way to check 
the quality of p, (/>-mesons as low-Q gluonometers 
would be to study color transparency in the pro- 
duction of these mesons at intermediate Q 2 , see 
discussion in section 3.7.1. 

Experiments at HERA have confirmed a num- 
ber of the QCD predictions [ J?J for the vector me- 
son production: (i) The rate of the increase with 
energy for p production for Q 2 > 20 GeV 2 and 



for J/tp production is xGn{x,Q 2 ^ j) 



oc W 



0.8 



for Ql ff ~ 4 GeV 2 . (Note that the soft physics 
expectation is ~ V^ 32 ). The models which take 
into account higher twist effects based on eq.(p"g|) 
[ |ll|, are able to to describe the HERA data at 
lower Q 2 . The inclusion of the quark exchange al- 
lows to reproduce the magnitude of the fixed tar- 
get data as well [ [l8| . (ii) absolute cross section of 
T production is reasonably reproduced [ [H]] . (iii) 
The extensive data of HI [ ^| on the energy and 
Q 2 dependence of J/ip production agree with pre- 
dictions of [ |ll[ and confirm need for large higher 

m 
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Figure 5. The dependence of average b and effective Q 2 of Q 2 for production of vector mesons [ [Tl| 



twist effects. The model of Ryskin [ £f) which 
neglects these effects leads to a factor ~ 4 — 5 
faster decrease of the cross section with Q 2 , (iv) 
the dominance of ol over (Tj- at Q 2 > m 2 , is 
confirmed (iv) The <f>/ p ratio has increased a fac- 
tor ~ 4 as compared to the value at Q ~ and 
reached the value close to the QCD prediction of 
1.2*(2/9) (where 2/9 is the SU(3) result), similar 
trend is observed for the J/tp/p ratio, (v) Con- 
vergence of the i-slopes of p meson production at 
large Q 2 and the J/ip slope (which weakly de- 
pends on Q 2 ) has been observed at HERA. The 
Q 2 dependence of B can be described semiquan- 
titatively as a result of the finite b in the vector 
meson wave function, see Figj6[ (vi) The data are 
consistent with a slower shrinkage of diffractive 
cone for these processes than for soft phenomena, 
see discussion in section 2.2. 
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Q [GeV 

Figure 6. Comparison the DESY data [ ^2| for 
the slope of the p meson production with the cal- 
culation of [ |TT| . 



2.6. BFKL dynamics - small dipole - small 
dipole scattering 

The cross section of scattering of two small 
dipoles was calculated in the leading log 1/x ap- 
proximation neglecting In Q 2 evolution by BFKL 

[§■ 

Calculation of next to leading log f /x effects 
were reported a year ago which demonstrated a 
large reduction of the rate of increase with energy, 
as compared to the leading log 1/x approxima- 
tion. To a large extent this reflected the smallness 
of the phase space - emission of one gluon in this 
approximation requires at least about two units of 



rapidity (at HERA this corresponds to emission 
of < 3 gluons). A number of attempts to perform 
resummation of NLO BFKL were reported re- 
cently. The current scenario is that a cx s - 25±0 05 
and that diffusion in the impact parameter space 
may be rather weak and not bring an onset of 
nonperturbative regime over a significant energy 
range 
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This scenario leads to a slower energy depen- 
dence of the cross section than eq.|f| for the scat- 
tering of a small dipole off a large hadron. Such a 
pattern seems unrealistic and reflects that in the 
HERA range the In Q 2 effects which are neglected 
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in the BFKL model are very important. Hence 
the quantitative comparison of the model with 
the data is hardly possible at present. It would 
be more promising to study the p t ordering of the 
produced hadrons which is qualitatively different 
in BFKL and DGLAP approximations. There- 
fore detection of particles over large rapidity in- 
terval is necessary - to emit three gluons between 
dipoles in BFKL kinematics Ay > 2 between the 
jets one needs a detector with acceptance of > 8 
units in rapidity. Promising direction of the study 
of the BFKL regime is 7* — 7* scattering at large 
Qi ~ Q2 ■ in this case contribution of the scatter- 
ing of two small dipoles is enhanced by a power 
of Q 2 as compared to the scattering of a small 
dipole on a large dipole [ [24) . 

2.7. Main focus of the future studies 

Since the most dramatic effects are expected to 
be happening in the processes directly coupled to 
gluons the major focus of the studies should be a 
measurement of the quantities more directly sen- 
sitive to the gluon dynamics and to the onset of 
hard (pQCD) regime. This includes measurement 
of <tl(x,Q 2 ), the absolute value, W,Q 2 ,t depen- 
dence of diffraction for the exclusive production 
both for c l and <jt (vector mesons, ttit ....), for 
charm, dijets, etc. 

The key issue is to establish nonuniversality of 
the vacuum exchange - the Q 2 and flavor depen- 
dence of the rate of the change of the i-slopes 
with energy, and of the energy dependence of the 
diffractive processes, and especially the difference 
between the hard diffraction induced by hard 7*- 
gluon and by 7*-quark scattering. 

Another direction of the research would be 
looking for manifestations of BKFL-type dynam- 
ics in the hard processes spanning maximal pos- 
sible range in rapidity. 

To achieve these objectives one would need to 
perform high precision DIS measurements both 
at lower energies down to W ~ WGeV and be- 
yond the energies of HERA - future ep colliders. 
Accurate measurements of the t dependence of 
diffractive processes requires improving the for- 
ward acceptance of the collider detectors. The ex- 
tension of the detector rapidity acceptance range 
much further into the proton fragmentation re- 



gion would be necessary to reach the rapidity 
ranges of Ay ~ 6 — 8 which would allow studies 
of a new pattern of pt ordering important in the 
kinematics where log j effects become essential. 

3. SMALL x PHYSICS IN eA SCATTER- 
ING 

3.1. High-Energy eA scattering - History 

Thirty years ago nuclei were suggested as a di- 
rect way to observe hadron properties of light, see 
review and references in [ . 

The striking effect was that in spite of the small 
value of ct 7 jv, the screening was observed for 7* 



-nucleus cross sections: 



Aa-., 



< 1. The reason 



is that a virtual photon transforms to a hadronic 
(quark- gluon) configuration \h) which propagates 
distances l coh = » in the DIS 

limit, which by far exceed the nucleus size of 
2Ra (large Ioffe distances). This key feature of 
space-time development of the high-energy DIS is 
a backbone of many analyses of small x DIS at 
HERA which we discussed in section 2. 

Hence there are two distinctive kinematics for 
DIS electron - nucleus scattering. At x > 0.2 
nuclear effects in the total cross section originate 
from the short-range nuclear structure, while at 
small x <C 0.1, l co h 3> 2Ra cross section is practi- 
cally not sensitive to details of the nuclear wave 
function. 

3.2. Structure functions of nuclei in the 
black disk limit 

In the limit l coh (x) > 2R A ( x < 10~ 2 for 
A=200) there is a deep connection between the 
7* A scattering and R e+e " (M 2 ) = a(e + e" -> 
hadrons) / 'a (e + e~~ — > fi + fi~) [ |J 



a_ r <T t ot( u h h 
3tt Af? 



A)R e 



(Ttot(l*A) 

~(M 2 )M 2 



{Q 2 + M 



2 \2 



dM 2 , (15) 



If one makes an assumption natural for the soft 
strong interaction dynamics that ALL hadronic 
configurations \h) in 7* at large enough v = 
E~* would interact with a large nucleus with 
o~tot("h" — A) ~ 2ttR 2 a a gross violation of 
the Bjorken scaling should take place for Q 2 = 
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const, v — > oo [ |j: 

w^ M2) h ? (16 » 

where xn = 75 — ^3— • 

Since the Bjorken scaling was considered a nat- 
ural behavior of the cross sections at large Q 2 this 
prediction of the gross violation of the Bjorken 
scaling was refereed by Bjorken as the Gribov 
paradox. Bjorken suggested the parton model so- 
lution of the paradox: the Aligned Jet Model - 
only a small fraction ~ 1/Q 2 of configurations 
in 7* which have small k t (large transverse size) 
could interact with hadron-like cross sections [ 
|2(| . The rest of configurations would not interact 
(parton model) or interact with in a color trans- 
parent way (pQCD of intermediate energies). 

In QCD, due to a fast increase with energy of 
the cross section of interaction of a "small color 
dipole" with a nucleon, the Gribov scenario be- 
comes ones again a viable possibility in the Q 2 
=const, qo — ► 00 limit for F2,Fl- In the case of 
gluons it leads to [ : 

x G A (x,Q 2 )/Q 2 = ^ln^. (17) 

Hence xG A {10' 4 ,Q 2 = 10) /A ~ which is 
about 6 for A=200 and should be compared 
to the current fits which give xGn(10~ 4 ,Q 2 = 
10) ~ 20. Thus eq.(|l7|) implies that at HERA a 
strong reduction of the gluon densities in heavy 
nuclei should be observed. The interaction in the 
gluon sector in the case of heavy nuclei would be 
close to the black body limit corresponding to a 
drastic violation of the DGLAP evolution. 

3.3. Summary of the HERA 95-96 eA 
study 

The study performed within the framework of 
the workshop "Future of HERA" [ ||| has identi- 
fied several fundamental questions which can be 
addressed within the electron-nucleus program: 
(i) amplification of the nonlinear effects expected 
in QCD at small x, (ii) the dynamics of high- 
energy interactions of small color singlet systems 
with nuclei, (iii) investigation of propagation of 
quarks and gluons through nuclear matter - en- 
ergy losses, pt broadening, extra hadron produc- 
tion, (iv) connection to the heavy ion physics. 



HERA would allow to use three complemen- 
tary tools: inclusive parton densities - quark and 
gluon nuclear shadowing, exclusive and inclusive 
hard diffraction, leading and central hadron pro- 
duction in inelastic eA collisions. Measurements 
with a set of nuclei with a luminosity 10pb~ 1 /A 
per nucleus would allow high precision studies of 
the dependence of various characteristics of DIS 
as a function of A at 100 times smaller x than in 
the fixed target experiments (which could study 
a small fractions of the key observables available 
for a collider). 

Frequent switching between different ions 
would lead to a strong reduction of the detec- 
tor related systematic errors for various nuclear 
ratios. Systematic errors due to the radiative cor- 
rections are also small. 

In particular, it would be possible to measure 
the difference the slopes d(FA /F 2 D )/d\nQ 2 and 
dF^/dlnQ 2 with accuracy ~ 5% (lpb _1 /A per 
nucleus) and measure directly Ga/Gn with ac- 
curacy of better than 10% (10-20pb _1 /A per nu- 
cleus). This would allow to observe 10% devia- 
tions from the DGLAP predictions. 

To achieve the desired luminosity it would be 
necessary to have luminosity decreasing approx- 
imately cx 4. With a reasonable capital invest- 
ment this condition can be almost satisfied for 
A ~ 16, while to reach this goal for heavy nuclei 
it appears necessary to perform cooling of heavy 
nuclei. In principle, the cooling may be easier 
for heavy ions than for protons. Its feasibility at 
HERA is currently is under investigation 

The development of the capabilities to accel- 
erate nuclei would make it possible also to per- 
form studies of the neutron parton densities using 
the deuteron beams and the proton tagging and 
hence to measure the nonsinglet parton densities. 
If the scheme for the accelerator of the polarized 
deuterons could be developed f\ this would allow 
to measure simultaneously the singlet and nons- 
inglet giN which would be very beneficial for the 
spin program. 



5 M.S. is indebted to F.Willeke and J. Maidment for dis- 
cussions of the accelerator related issues. 
6 V.Skrinsky, private communication 
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3.4. Infinite momentum frame approaches 
to nuclear shadowing 

There exist two complementary approaches to 
descriptions of the small x physics: one is based 
on the use of the parton frame (fast nucleus) 
another on the use of the nucleus rest frame. 

The fast frame descriptions are convenient for 
the pQCD analyses, modeling of nonlinear effects. 
They usually assume that there is no shadowing 
in the leading twist at Qq [ [2£ 



30, 31 P2 33 



The rest frame descriptions allows to visual- 
ize connection to the soft physics and to the ep 
diffraction [ [§, |§ |§ ^ || |3§ g || £| g . 
So far the effects of the higher twist effects were 
not treated explicitly in these approaches. 

In the fast frame nucleon looks as a "MA- 
TRESHKA ": the fast partons in nucleons are con- 
tracted to a pancake of radius Ra and a small 
longitudinal size z — 2tm^ m - while the small x 

partons form pancakes of times larger size 
{(x v ) ~ 0.2 is average x of the valence quarks). A 
nucleus looks as a collection of MATRESHKA's 
with valence quarks not overlapping and clouds 
of partons of different nucleons with x <C (x v ) -In- 
completely overlapping. 

One can consider small x parton density per 
unit transverse area: ^pqRaxGn(x,Q 2 ), where 
Po r~j 0.16/to~ 3 is the average nuclear density. 
At central impact parameters parton density/unit 
area is 1.5 times larger: 2pqRaxGn(x,Q 2 ). 
Hence one finds a large enhancement of the 
parton densities/area in nuclei as compared to 
nucleons. For the central impact parameters: 



G A /irR 2 A 



Ai 



Ga 



G 



1.4=200 



11 AGn 1 



G N /nr% ~ AGn 

One can study interaction of overlapping tubes 
of partons in the perturbative strongly interact- 
ing phase In the simplest model of Mueller and 
Qiu [ [2{| nonlinear effects are described by the fan 
diagrams. The additional nonlinear term contri- 
bution SGa{x, Q 2 ) in the evolution equation for 
G A (x,Q 2 ): 

d SxG A (x,Q 2 ) _ 



81 A 1 / 3 
16 Q 2 r 2 



o? s (Q 2 ) 



dQ 2 A 

1 —[uG N {u,Q 2 )] 2 . (18) 



O 
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to 
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For A ~ 200 this correction is three times larger 



Figure 7. The scaling violation for the nuclear 
gluon density for A=200 calculated in the leading 
twist and with inclusion of the Muller, Qiu term. 



than for the proton and would strongly reduce 
a t Q2 = 4 GeY 2 x „ 1Q -3 _ 10 -4 )A = 

200 as compared to the DGLAP predictions 
(Fig.7). 

Note however that if this correction is large 
it certainly cannot be trusted since other effects 
should become important. Current models of 
nonlinear effects try to account for multinucleon 
effects, etc. They involve many simplifying as- 
sumptions. The overall trend is an expectation of 
large gluon shadowing even if starting from the 
assumption that at a low normalization point the 
leading twist shadowing is small. 

Though this approach is seemingly qualita- 
tively different from the rest frame approach 
which we started with and which we will discuss 
in the next subsection these descriptions are in 
fact dual. 

When viewed in the rest frame deep inelas- 
tic scattering corresponds to the scattering off 
a small dipole off the nucleus. The dipole re- 
solves strong gluon fields at the scale Q ~ t and 
propagates much longer distances 2Ra versus 2tn 
through the gluon fields leading to amplification 
of nonlinear effects. Therefore large nonlinear ef- 
fects are expected based on the space-time pic- 
ture in both reference frames both in the struc- 
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ture of the final states and in a very large gluon 
shadowing which should evolve with Q 2 in a way 
substantially different from the DGLAP expecta- 
tions. 

3.5. The leading twist shadowing for par- 
ton densities and hard diffraction in 
electron-proton scattering 

There exists a deep connection between high- 
energy diffraction and phenomenon of nuclear 
shadowing [ || . Qualitatively it is due to a pos- 
sibility of small momentum transfer —t m i n = 



Xjpm 



N 



where xp = Xbj(l 



M dif/Q 2 )- IfV=*< 



"average nucleon momentum in D(A)" the am- 
plitudes of diffractive scattering off proton and 
off neutron interfere. The corresponding double 
scattering diagram for the j*D scattering leads to 
a diffractive final state. The Abramovskii, Gri- 
bov, 
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then implies that 



Kancheli theorem 
this interference for a deuteron or a light nu- 
cleus: (i) increases the diffractive total cross sec- 
tion by Audi f = 02, (ii) decreases total cross sec- 
tion by Aatot = —<?2, (hi) decreases cross sec- 
tion of inelastic interactions with a single nucleon 
by Aa a i n gi e = — 4<72> (iv) results in simultane- 
ous inelastic interactions with two nucleons with 
o 'double = 2<72- The sum of partial contributions 
gives the contribution of the diffraction to the to- 
tal cross section: Aatot(eN) — a to t{eA) = Aer tot = 
Audif + A 

^single T" ^double- 

Thus there exists a direct connection between 
Cdi//(eiV) at t ~ and A<j t ot(eA) for light nuclei. 
Neglecting small corrections due to the real part 
of the diffractive amplitude one finds: 

da diJf (ep) 
ed _ dt l*=0 



dt 

8ttR 2 d 



(19) 



where Rd is the deuteron radius. 

As we discussed in section 2 the QCD factor- 
ization theorem is valid for DIS diffraction: the 
Q 2 evolution of diffractive structure functions at 
fixed xp,pt is described by DGLAP. By compar- 
ing the QCD diagrams for hard diffraction and 
for nuclear shadowing due to scattering off two 
nucleons (Fig. 9) one can prove [ that in the 
low thickness limit the leading twist nuclear shad- 
owing is unambiguously expressed through the 
diffractive parton densities fP(^f-,Q 2 ,xp,t) of 



ep scattering: 

fj/A(x,Q 2 )/A = f j/N (x,Q 2 ) 

^ r r°° r°° r x o 

— / d 2 b / dz\ I dz2 I dxp ■ 

2 J J —GO JZl J X 

■fj/N (/3,Q 2 ,£iP)*)| fe 2 =0 PA(b, zi)p A (b, z 2 ) ■ 

cos(xpm N (z 1 — z 2 )), (20) 

where Q 2 ), fj/ N (x, Q 2 ), are inclusive par- 

ton densities; Pa{t) is the nucleon density in the 
nucleus. At not too small x one should add a term 
related to the longitudinal distances comparable 
to internucleon distances in nuclei. This addi- 
tional term can be evaluated using information 
on enhancement of the gluon and valence quark 
parton densities at x ~ 0.1 in the normalization 
point as an input. It slightly diminishes nuclear 
shadowing at higher Q 2 via the Q 2 evolution. 
Thus for light nuclei at small enough x the shad- 
owing effect for 1 — f 3 /A(x,Q 2 )/Afj/ N (x,Q 2 ) cx 



crf // (x,Q=2) 
7« ff (x,Q=2) 



ai ff (x,Q 2 ) = 16nda di f/dt lt=0 /a tot . 

As we discussed before, the analyses of 
the HERA hard diffractive data indicate that 
a 9 eff (x < 10~ 3 ,Q = 2) - 50m&, 
3.5. 

If one follows indications from HERA that in 
the Q 2 > 4 GeV 2 range the leading twist already 
dominates in diffraction (this is well established 
for quarks, but may not be a good approximation 
for the gluons (see discussion in section 2.2) one 
predicts: (i) a significant shadowing in the quark 
channel with small theoretical uncertainty, (ii) a 
much larger (a factor of 3 for light nuclei) shadow- 
ing in the gluon channel subject to the question 
of the higher twist effects (Fig.^). 

Note than in many models of nuclear shadow- 
ing which did not incorporate information about 
the gluon induced diffraction it was expected that 
the shadowing in the gluon channel would be 
smaller than in the quark channel. 

The Gribov theory agrees well with the high 
precision NMC data if one uses the HERA diffrac- 
tive data (this involves a certain extrapolation of 
the HERA data to smaller W) [ g| . A study of in- 
clusive parton densities at the HERA x, Q 2 range 
would (i) establish shadowing in the quark chan- 
nel, and check connection between diffraction and 




Figure 9. The dependence of G A /AG N and F 2 a/AF 2 n on x for Q=2 ( dashed), 5 (dotted), 10 GeV ( 
solid) curves calculated in the quasieikonal model using diffractive parton densities of [ H. 
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Figure 10. Suppression of the jet production in 
AA collisions due to gluon shadowing at y = 
calculated in the quasicikonal and fluctuation 
models (solid and dashed curves). 



shadowing starting from nonperturbative domain 
of Q 2 ~ where shadowing strongly depends on 
Q 2 , (ii) find out whether interactions in the gluon 
channel are superstrong (~50 mb) at HERA for 
intermediate Q ~ 2 — 3 GeV, (iii) study coherent 
interactions of 7* with three, four, ... nuclcons in 
the scattering off heavy nuclei. Such interactions 
are sensitive to the fluctuations of the interac- 
tion strength in the quark and gluon channel and 
check the role of weakly interacting components - 
this information cannot be obtained from ep scat- 
tering. For A=200 fluctuations may reduce the 
gluon shadowing at x = 10~ 3 — 10~ 4 by a factor 
up to ~ 1.5 [ |2?j]. Hence the study of shadow- 
ing for A > 40 would provide a complementary 
information to the case of light nuclei; (iv) verify 
connection of diffraction and shadowing via inves- 
tigation of the Q 2 dependence of shadowing be- 
tween Q 2 = and 4 GeV 2 which should be strong 
due to a strong Q 2 dependence of the diffraction 
contribution to cr 7 *Ar. 

3.6. Implications for the heavy ion colli- 
sions at LHC 

The production of minijets is often considered 
as an effective mechanism of producing high den- 
sities in the head on heavy ion collisions. For 
the central rapidities minijets are produced due 




10 " J 10 " 3 10 ' 2 



Figure 11. Color opacity effect for the ratio of 
the coherent production of J/tp from carbon(lead) 
and a nucleon normalized to the value of this 
ratio at x = 0.02 (where color transparency is 
expected to be valid) calculated in the leading 
twist approximation with (dashed) and without ( 
solid) account of the fluctuations of the interac- 
tion strength. 



to collisions of partons with £ 7et = For 
heavy ion collisions at LHC y^ijviv ~ 4 TeV and 
the gluon-gluon collisions are responsible for pro- 
duction of most of the minijets. Therefore the 
gluon nuclear shadowing would lead to a reduc- 
tion of the rate of the jet production due to the 
leading twist mechanism by a large factor up to 
Pt ~ WGeV/c, see Figjl^ where we give results 
of calculation in the quasieikonal and fluctuation 
models of shadowing. The nuclear gluon shad- 
owing leads to a similar very strong reduction of 
the heavy onium production in pA and AA colli- 
sions at LHC energies for y c .m. 

— and small p t . 
However nuclear shadowing leads to other contri- 
butions to the jet production not included in the 
factorization models - fragmentation of the spec- 
tator partons belonging to the dipoles, decay of 
the ladders involved in the multiple interactions 
with the dipoles. Hence the overall effect could be 
a smaller suppression and additional emission of 
soft particles over a large rapidity range. Exper- 
imental studies of eA collisions would be of great 
help in quantifying these effects. 

3.7. Study of diffraction with nuclei 
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3.7.1. Exclusive channels: M = p, J/ip, irir, ... 

The QCD factorization theorem for exclusive 
processes [ |l6| leads to [ |?J 

w(rA^MA)\ t=Q 



w(rN^MN)\ t=Q 



r F%(x,Q) 
Fk(x,Q) 



G\(x,Q) 
G N (x, Q) 



(21) 



Hence in the region x > 0.01 the color trans- 
parency regime of weak absorption is expected. 
Establishing down to what Q 2 this regime would 
hold would allow to determine down to what vir- 
tualities production of vector mesons could be 
used as a low Q 2 gluonometer. 

At smaller x the gluon shadowing would lead 
to a dramatic decrease of the yield of vector me- 
son production at small x (yields are still large 
to measure the cross section accurately) - the 
color opacity phenomenon - strong absorption of 
small color dipoles by nuclei at ultra-high ener- 
gies. This phenomenon is a qualitative departure 
from the experience of QED where cross section 
of interaction of small electric dipoles practically 
does not depend on the incident energy. The pre- 
dicted A-dependence for coherent J / ip production 
is about the same as for photoproduction of p- 
mesons!? 

3.7.2. Inclusive and semiinclusive chan- 
nels: 

Prime objectives of the study of the inclusive 
diffraction is to obtain a direct answer to the 
question how "black" are diffractive interactions. 
We expect an increase of the rapidity gap prob- 
ability with A. For the heavy nuclei the fraction 
of events where a nucleus remains intact would 
reach ~ 30% for the quark induced hard diffrac- 
tion and nearly 50% for the gluon induced diffrac- 
tion [ fil], [27]]. Also, the rapidity gap probability 
would depend on p t of the jet is a rather pecu- 
liar way - contrary to naive expectations it would 
increase between the aligned jet kinematics and 
the pt range where gluons give dominant contri- 
bution. 

Overall the study of inelastic diffraction would: 
provide a direct test of the strengths of interac- 
tion in the quark and gluon sector - identify the 



channels for which interaction is close to black 
for heavy nuclei. This information will be com- 
plementary to the study of the A dependence of 
Ga/Gn and allow to clarify the origin of gluon 
component of the "Pomeron'. 



3.8. Nondiffractive hadron production 

Numerous data on hadron-nucleus scatter- 
ing at fixed target energies indicate that the 
multiplicities of the leading hadrons Na(z) = 
1_ da ( z ) decrease with increase of A. 

crtot(aA) dz 

Here z is the light-cone fraction of the projec- 
tile " a" momentum carried by the hadron "h" . 
On the contrary, the QCD factorization theorem 
for the inclusive hadron production in DIS im- 
plies that in the case of electron-nucleus scat- 
tering no such dependence should be present in 
the DIS case. This indicates that there should 
be a interesting transition from the soft physics 
dominating in the interactions of real photons 
with nuclei to the hard physics in the inclusive 
hadron production in the DIS kinematics. It 
would be manifested in the disappearance of the 
A-dependence of the leading spectra at large z: 
N A (z,Q 2 ) = N N (z,Q 2 ), for z > 0.2, Q 2 > few 
GeV 2 . At the same time the transverse spectra 
are expected to be broadened for nuclei (Ap 2 oc 
A x / Z ) due to the QCD analog the Landau-Migdal- 
Pomeranchuk effect, see discussion and references 
m[@. 

At small x a new interesting phenomenon 
should emerge for smaller z due to the pres- 
ence of diffraction and nuclear shadowing. In- 
deed, the diffraction originates from the presence 
in the wave function of 7* of partons with rela- 
tively small virtualities which screen the color of 
the leading parton(partons) with large virtuality 
and can rescatter elastically from a target (sev- 
eral target nucleons in the case of nuclear tar- 
get). Inelastic interactions of these soft partons 
with several nucleons should lead to a plenty of 
new revealing phenomena in small x DIS eA scat- 
tering, which resemble hadron-nucleus scattering 
but with a shift in rapidity from y max (current) 
related to the average rapidities of these soft par- 
tons. This shift can be expressed through the 
average masses of the hadron states produced in 
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Figure 12. A sketch of the A-dependence of the 
inclusive hadron production in DIS as a function 
of rapidity. 



the diffraction: 

Vsoft partons ~ ~ ]n((M^f) /[J 2 ), (22) 

where \i ~ lGeV is the soft scale. Partons 
with these rapidities will interact in multiple col- 
lisions and loose their energy leading to a dip in 
the ratio n A (y) = A A (y)/A p (y)(Fig.l2). At the 
same time these multiple interactions should gen- 
erate a larger multiplicities at smaller rapidities. 
Application of the AGK rules indicates that for 
V < Vsoft partons- A, where A = 2-3 the hadron 
multiplicity in the case of nuclei will be enhanced 
by the factor: rj A (y) = AF 2p {x,Q 2 ) / 'F 2A {x 1 Q 2 ). 
At the rapidities close to the nuclear rapidities a 
further increase of T)a{u) is possible due to forma- 
tion of hadrons inside the nucleus. 

One also expects a number of phenomena due 
to long range correlations in rapidity. This in- 
cludes: (a) local fluctuations of multiplicity in the 
central rapidity region, e.g. the observation of a 
broader distribution of the number of particles 
per unit rapidity, due to fluctuations of the num- 
ber of wounded nucleons [ [ll| . These fluctuations 
should be larger for the hard processes induced by 
gluons, for example the direct photon production 
of two high p t dijets. (b) Correlation of the cen- 
tral multiplicity with the multiplicity of neutrons 



in the forward neutron detector, etc. 

To summarize, the small x eA physics is one of 
the last unexplored frontiers of QCD. A plenty of 
new nonlinear QCD phenomena are to be discov- 
ered in these studies: the gluon nuclear shadow- 
ing, the high-energy color transparency, the per- 
turbative color opacity,.. These are just expected 
discoveries - what about unexpected ones? Studies 
of eA collisions would provide decisive tests of the 
interface of pQCD and soft physics, allow to reach 
understanding of the space-time picture of the 
high-energy strong interactions. Important appli- 
cations of these studies include the measurement 
of the nonsinglet parton densities in eD scatter- 
ing, as well as numerous connections to the RHIC 
and LHC heavy ion physics. 

HERA has produced exciting glimpses of new 
small x strong interaction gluon dynamics. The 
aim now is to produce a definitive proof and ob- 
serve high density hard QCD. 
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